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Various experiments evidence spatial heterogeneities in sheared yield stress fluids. To account 
for heterogeneities in the velocity gradient direction, we use a simple model corresponding to a 
non-monotonous local flow curve and study a simple shear geometry. Different types of boundary 
conditions are considered. Under controlled macroscopic shear stress E, we find homogeneous flow in 
the bulk and a hysteretic macroscopic stress - shear rate curve. Under controlled macroscopic shear 
rate V, shear banding is predicted within a range of values of T. For small shear rates, stick slip can 
also be observed. These qualitative behaviours are robust to changes in the boundary conditions. 

PACS numbers: 83.60La,83.10Tv,83.10Gr 



T3 

C 

O 

> 
O 

(N 

o 

(N 
O 
-i— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

Complex materials, such as concentrated colloidal sus- 
pensions, gels or foams, often exhibit exotic rhcological 
behaviour at low shear @, § § § | § f| § g, ^ [H), 
very different from the Newtonian linear response char- 
acterised by a stress £ linearly proportional to the shear 
rate T (in a scalar language). A term often coined for 
some of these behaviours is that of "yield stress fluids" . 
Although often loosely used, it most of the time describes 
fluids that have one or both of the following properties: 
when one looks at the steady state flow curve £(F) at 
low shear, i) the stress E tends to a finite value £y when 
the shear rate goes to zero, ii) no flowing steady state 
solution exists for stresses smaller than a value Eo- Con- 
sequently at low stresses, the system does not really flow, 
and these "pastes" often display physical aging in a way 
rather similar to structural glasses. 

A traditional way to apprehend this macroscopic rhe- 
ological non-linear behaviour is to suppose the existence 
of a local relation between the local stress a and the lo- 
cal shear rate 7 of one of the forms I to IV on figure [j]. 
Classically again, I to III being monotonic can lead to 
a homogeneous flow so that the macroscopic response is 
identical to the local flow curve : while II and III cor- 
respond to genuine yield stress fluids satisfying both i) 
and ii) above (with Eo = Ey = (7(7 = 0) ), I yields a 
power law fluid behaviour E = AT a which can be exper- 
imentally very difficult to tell apart from the yield stress 
behaviours II and III if a is small. It is on the other 
hand classically argued that the decreasing branch of IV 
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FIG. 1: Schematic representation of possible local steady 
state curve cr(j) generating a macroscopic yield-stress fluid 
behaviour. 



(dYj/dT < 0) can not lead to the corresponding homoge- 
neous situation, as such a flow is unstable with respect 
to the formation of "shear bands" or domains flowing at 
different local shear rates. This statement is usually sub- 
stantiated by a simple argument invoking inertial terms 
(see, for instance fL2|). This type of banding is referred to 
as 'gradient banding' by opposition to 'vorticity banding' 
(coexistence of domains flowing at different local shear 
stress and common shear rate). 

The question then arises of the structure of the 
heterogeneous flow pattern observed in these conditions 
(curve IV, gradient banding), of its stability (i.e. are 
there always steady-state solutions?), and of the result- 
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ing macroscopic flow curve £(r) for the steady-state 
solutions. Investigating this problem in generic terms 
is very difficult as many underlying dynamical laws 
involving an unlimited number of "hidden" variables can 
generate local equations such as IV. In this paper, to 
make progress, we consequently make a specific choice 
for this underlying dynamics, that we justify now. 

Formally, taking the vertical 7 = axis as part of 
the local flow curve, the problem is very similar to 
the analysis of flow-induced or flow-assisted "phase" 
transitions for systems with a local non-monotonic 
relationship (7(7). For some systems, the change in 
rheological properties can be ascribed to a structural 
phase transition occurring under flow : Ramos et al. 
with lyotropic hexagonal phases Q], Diat et al. with 
onions [|] and Britton et al. and Berret et al. with 
wormlike micelles |j| observe a phase transition, 
and the coexistence of two well defined phases. Such 
behaviours are quite reminiscent of equilibrium phase 
transitions, although no general principle applies a 
priori in such out of equilibrium situations ]lq ]. The 
identification of structural differences between the flow- 
ing phases then provides a guide as to what are the most 
likely relevant local variables beyond the rheological 
ones: a (non-conserved) nematic order parameter, the 
(conserved) concentration of a constituent, etc... Related 
theoretical descriptions have thus tried to describe the 
coupling between this order parameter and flow in 
order to obtain a selection criterion for the flow pattern 
and the resulting macroscopic rheological behaviour 

However for many yield stress fluids there is no ob- 
vious structural difference between the coexisting flow- 
ing phases: experimental reports often describe "shear 
bands" of a flowing "phase" in a stationary "phase" with 
no clear structural or textural difference between the two 
parts of the system. Such "soft glassy materials" g in- 
clude gels, pasty systems, concentrated glassy suspen- 
sions, . . . One may cite the work by Raynaud et al. [Q, 
where banding is observed in a suspension of bentonitc 
and the interface between the bands is in the velocity- 
vorticity plane (gradient banding) : the system is sheared 
within a band at the inner moving cylinder, while the rest 
of the sample remains un-sheared. Laponite suspensions 
also undergo banding for a certain range of imposed shear 
rates, range for which the shear stress is almost constant 
|jTof . Noticeably, for a range a lower shear rates 'stick 
slip' is observed in this system : a layer reversibly frac- 
tures and re-heals, showing stationary oscillations in the 
measured stress. Eventually, one may note that inho- 
mogeneous velocity profiles are also observed in sheared 
foams Jf| , and granular pastes . 

For these systems with no clearly identified structural 
variable that would distinguish coexisting phases, there 
is much less guidance for building models. Still, banding 
is likely to originate from the same basic physical mecha- 



nism as for the 'phase separating' fluids (20) . Models have 
therefore been derived along similar lines, which differ by 
the choice of the hidden variable (equivalent to the or- 
der parameter) and of its dynamics (20[ ^2). In addition 
or alternatively, gradient terms were included in the non 
monotonic intrinsic law (7(7) [fl2| , plf . 

In the present paper we follow the route initiated by 
Derec et al. p3[ , that has been shown to be a simple 
though useful guide to understand the interplay between 
non-linear rheology and aging properties for yield-stress 
fluids characterised by local flow curves of types I to III 
[2~|] . In the present work we adapt it to describe local 
flow curves of type IV, and to take into account spatial 
heterogeneities. 

The essence of the model of [^3) is to describe the local 
state of the system by a single scalar variable that is for 
convenience taken to be the local relaxation frequency of 
the stress (called "fluidity"), as expressed by a Maxwell 
equation: dt& = —a(t)o~ + 7 (stress has been rescaled 
so that the elastic modulus is 1). Then the evolution of 
the fluidity a(t) results from the competition between: (i) 
aging - relaxation rates in the system decrease during the 
spontaneous evolution of the pasty system, this evolution 
being itself slower and slower as time passes; and (ii) 
flow-induced rejuvenation -, flow induces rearrangements 
in the system, which in turn increase the "fluidity" of 
the system as a feed back mechanism || ^4|, |2{| . This 
competition results in an equation for a that supplements 
the Maxwell one: dta = —/(a) + g(a, er, 7), with / and g 
positive functions. 

Here, to address the question at hand, we focus on 
choices for / and g leading to local flow curves (local 
steady-states) with a decreasing branch at low values of 
7, and we attempt to describe spatial heterogeneities, 
through a simple diffusion term DAa in the equation 
describing the evolution of a. 

We then study the flow patterns in a simple shear ge- 
ometry for such a system, when it is driven at cither 
imposed macroscopic shear stress £ or imposed macro- 
scopic shear rate T, and deduce the resulting macroscopic 
rheological behaviour. We will on the way provide ana- 
lytical and numerical arguments, using when necessary 
explicit choices for / and g, and pinpoint the role of the 
boundary conditions. In particular we will show that the 
model accounts for the following features : 

• when £ is imposed, the flow curve jumps in a hys- 
teretic fashion between two branches: one that cor- 
responds to no flow but for a wall layer in the vicin- 
ity of the walls; the other corresponds to a fully 
fluidised situation. 

• when the global shear rate T is imposed, shear 
banding can occur, as well as sometimes a stick-slip 
like oscillating behaviour at small shear rate that 
corresponds to a localised oscillation of the fluidity. 

These features will be analysed in the context of experi- 
mental results. 
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The paper is organised as follows. We recall in Section 
H the basic ingredients of the model in Ref. ||| and 
propose a generalisation to account for inhomogeneities 
in the system. The homogeneous case is detailed in sec- 
tion III and the stability of the homogeneous solutions 
is discussed. Heterogeneous solutions are considered in 
Sections IV and[v|. In section IV, the controlled stress sit- 
uation is considered, while in section |v|, the global shear 
rate is imposed. In section VI, choices made in building 
the model and the influence of the boundary conditions 
are discussed. 
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FIG. 2: Schematic drawing of the local steady state flow curve 
f(7), °(a),j(a)- 



II. MODEL AND GEOMETRY 

In this section we define more clearly the model and 
the geometry we are going to study, following the physical 



arguments given in the introduction. In subsection II A, 
we s pecify further the functions / and g, in subsection 
II B we describe the geometry considered and the corre- 
sponding description of heterogeneities in the equations, 
in subsection II C we define the boundary conditions for 



both velocity and fluidity, and we summarise in subsec- 
tion [ID by recalling the resulting set of equations. 



A. Accounting for a non-monotonous local flow 
curve 

From the introduction, we start with the model of [B3j , 
devised for homogeneous situations: 



dt& = —aa + 7 

d t a= -f(a)+ g(a,(7,j) 



(1) 



For sake of explicitness, and at loss of generality, it will 
often be useful to restrict the choices for the functions / 
and g appearing in equation (1). This is in particular true 
for the flow induced rejuvenation term g(a, a, 7) which 
can take many forms. 

Following again Derec et al. p3]| , we take forms cor- 
responding to low values of the stress, shear rate and 
fluidity, so that power laws may express a kind of expan- 
sion of more general functions. Also a being a positive 
scalar, the direction of shear is of no influence so one ex- 
pects g(a, a, 7) = g{a, — cr, —7). A simple solution is then 
g{a,a,j) = h(a)aj. 

Then following the introduction, we want the steady- 
state flow curve <r(j) to have a decreasing branch for 
small but non-zero values of 7 (as sketched in figure 2). 
In that case, the steady-state solutions of the system (1), 
consist of : the vertical axis ( a — 0, 7 = 0, a arbi- 
trary), and solutions given by cr(a) = (/(a)/a/i(a)) 1 / 2 
when expressed in the plane s, a. The function 0(7) is a 
monovalued increasing function, whereas c(a) is multi- 
valued, decreasing for small values of a, then increasing 
for larger values. The requirements are thus that that 
f(a)/ah(a) has a negative derivative for < a < a m , 
a positive one for a m < a. The minimum thus reached 



for a m corresponds to the minimal value of the stress 
o"m = (j "(flm) / 1 0'mh{a m )) 1 / 2 . Eventually, as for curve IV 
in Figure 1, a tends towards a finite value when 7 goes 
to zero, which we call &m- 

When numerical calculation will prove necessary to 
analyse the evolution of the system, we make even more 
explicit choices regarding / and h, that need of course 
satisfy the above requirements. We use the simple choice: 
/(ft) = ra 2 + a 4 , h(a) — r'a + a 2 with r and r' arbitrary 
number that we vary to check the robustness of our con- 
clusions. That /(ft) scales as a 2 for small values of a 
gives a spontaneous relaxation of a as the inverse of the 
elapsed time, a situation often referred to as "full aging" 
(see [p3|). 



Incorporating heterogeneities: geometry and 
model 



In principle, a general description of heterogeneous 
flows implies complex algebra describing the evolution of 
scalar and tensorial quantities through convected deriva- 
tives, and writing various conservation equations (for 
mass, momentum, etc ..). We focus here on the simplest 
flow geometry (simple shear), allowing only variations 
along one direction, perpendicular to flow and vorticity. 

More precisely, we consider our complex fluid to be 
bounded by two parallel plates located at z — and 
z = H (see e.g. Figure 3). The top plate moves along 
x at a velocity V, while the bottom one is immobile, so 
that the macroscopic shear rate is T — V/H. In this 
situation we consider laminar flows along x, where the 
local variables a = a xz , 7, and a, depend only on z (i.e. 
not on x and y). 

Neglecting acceleration (inertial) terms is a common 
assumption for such (highly viscous) fluids on the ground 
of a low Reynolds number. This in general makes the 
stress divergence free, but in the present geometry it 
results in the stress being spatially constant so that 
it is no more a local variable. We will in this paper 
make this assumption of fast momentum diffusion, while 
keeping in mind that by doing so we formally suppress 
a "mechanical instability" that destabilises into bands 
any homogeneous flow solution corresponding to a de- 
creasing part of the local flow curve da/dj < (see 
e.g.p^, 20, |27|]). In this framework, the local stress de- 
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appearing from the competition of /(a) and Dd zz a. If 
- 1 as in the example at the end of the previous subsection 

-tj' f(a) ~ a 2 for small a, then this length scale depends on 

the actual value of a : (D/a) 1 / 2 , indicating algebraic tails 
~j for spatial structures of low fluidity. 



FIG. 3: Geometry: simple shear along x, variations along z 
only. The flow can be heterogeneous, the example depicted 
here correspond to a low fluidity in the bulk and wall layers 
of higher fluidity close to the wall. 



pends solely on time, and is identical to the macroscopic 
stress a(z,t) = <r(t) = £(i). 

Now in order to discuss spatial heterogeneities, we need 
to account for the effect on the dynamics of a heteroge- 
neous value of the fluidity. Again for the sake of simplic- 
ity, we model this as a diffusion process in the evolution 
of the fluidity and correspondingly add a DV 2 a diffusion 
term to the second equation of (1). This choice corre- 
sponds to a picture in which the system is fluidised in 
a localised region, this triggers spontaneous local relax- 
ations which themselves tend to fluidise the neighbour- 
hood of the initial region. Accordingly we take the "dif- 
fusive coefficient" to be a positive constant. A negative 
value would generate local, small wavelength instabili- 
ties that destabilise homogeneous flow even at high shear 
rates. As we focus here on the instability generated by 
the non monotonous branch of the local flow curve, and 
in line with experimental observations |Q, ^) of homoge- 
neous flow at high shear, we discard this possibility here. 
In addition, we do not consider other terms formally al- 
lowed by symmetry such as (Va) 2 . 

The system is therefore described by the coupled equa- 
tions: 



C. Boundary conditions 

The equation for the fluidity and velocity have to be 
supplemented by boundary conditions at the confining 
walls at z — and z = H. 

We here exclude any direct wall slippage and assume 
the plates to be sufficiently rough for the plate velocity 
to be equal to the limit velocity in the fluid approaching 
the plate (no slip velocity). This is equivalent to: 




dzj(z,t) = V(t)=t(t)H (3) 



It is important to point out that this hydrodynamic con- 
dition is not sufficient : a boundary condition for the 
fluidity equation is also required. How the walls affect 
the density, order, relaxation processes of the fluid in 
their vicinity is obviously a complex problem with many 
aspects to it, even for simple liquids p9|. A smooth flat 
wall can enhance ordering, a rough or chemically disor- 
dered one can on the contrary locally disorganise the fluid 
(relative to what happens in the bulk) . The consequences 
of these on the local relaxation rates being far from triv- 
ial p(J , we choose here to investigate two limiting simple 
cases that hopefully bracket the true behaviour : (i) a 
fixed fluidity at the wall 

a(z = 0) = a(z = H) = a w (4) 

and (ii) a zero gradient condition 

d z a\ z =o lH = 0. (5) 



f d t a{t) = -a(z,t)a(t) + j{z,t) 

) d t a(z,t) = -f(a(z,t)) + h(a(z,t))a(t)j(z,t)+ (2) 
{ Dd 2 zz a(z,t) 

One may note that in its present form, the equation for 
the fluidity is close to the equation for the granular tem- 
perature obtained within granular hydrodynamics of in- 
elastic beads |^8j : the spatial diffusion term accounts in 
the latter case for " heat" diffusion (according to Fourier 
law) and the " aging" term corresponds to the energy dis- 
sipation term (due to the inelasticity of collisions) , which 
drives the granular material towards a state of zero tem- 
perature. Eventually, the granular temperature is indeed 
a measure of the local fluidity of the granular material. 

An important comment is also that the diffusion coef- 
ficient D that we took constant here could in principle 
be a function of a. This would affect the length scale 



D. Summary 

We summarise here the equations that we will use in 
the remainder of the paper. 

' 9 t E = -aE + 7 
d t a = -f(a) + h{a)^ + Dd 2 zz a 

' J" dzj(z,t) = V = TH () 
t a\ z =a,H = a w or d z a\ z=0 ^ H = 

where S(t), T(t) are global variables whereas a(z,t) 
and j(z, t) are local ones. When explicit choices will be 
necessary to study the corresponding dynamics, we will 
use /(a) = ra 2 + a 4 and h(a) = r'a + a 2 . 
Let us recall that by neglecting the inertial terms, we 
formally suppress the "mechanical instability" that 
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destabilises the decreasing branch da/d-f < 0. 

Numerical method - For various cases, it is diffi- 
cult/impossible to get analytical results and numerical 
calculation becomes necessary. The evolution equations 
for the stress and the fluidity, Eq. (||), have been solved 
numerically using a pseudo-spectral method, together 
with a Runge-Kutta integrator This technique 

integrates the diffusion part of the time equation for 
the fluidity in Fourier space, while the non linear part 
is solved in real space. Special care has been taken to 
avoid aliasing of the solutions. A standard Fast Fourier 
Transform (FFT) algorithm is used for the forward 
and backward transformations between real and Fourier 
space. 

We wish to use equations (6) to deduce the macro- 
scopic behaviour of the system, to analyse the conditions 
for the existence of a steady-state, and to describe the 
features of the flow. Note that the boundaries appear 
explicitely, as potential heterogeneities, so that most of 
the time the structure close to the boundaries will be 
different than in the bulk, allowing the walls to nucleate 
new "phases" before the bulk. Actually the absence of 
noise terms in these equations suppresses the equivalent 
of homogeneous nucleation in thermodynamics problems. 
We briefly comment upon the addition of noise when it 
plays a significant role. In the following three sections, 
we first revisit briefly the formal case of a homogeneous 
behaviour (section III), we then analyse the simpler case 
of the response to a fixed stress E (section IV), before 
turning to the more complex case of a fixed global shear 
rate T (section V). 



III. HOMOGENEOUS BEHAVIOUR 

We first analyse briefly the homogeneous case, where 
we impose that the fluidity a(t) is independent of z, so 
that E and T are similar to a and 7. Although slightly re- 
dundant with previous discussions of the flow curve, this 
analysis will highlight the potential differences between 
the stability analyses at fixed stress and fixed shear rate. 
These differences naturally result from the fact that E 
and f are treated on a different footing in equation (0). 

As previously stated, steady-state solutions consist of 
three branches: 



branch "0" : r = 0, a = 0, E arbitrary. This 
corresponds to a frozen, non-flowing system. 

branch "I": a < a m , T < j m = a rn a m , and E a 
decreasing function of both a and T, 

branch "2": a > a m , V > 7 m , and E an increasing 
function of both a and T, from a rn up to 00. 

The three solutions (or branches) of fluidity arc schemat- 
ically plotted in the (E,a) plane on figure (the E(r) 
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FIG. 4: Schematic description of the evolution of a at fixed 
stress (horizontal arrows). Stable steady-state branches are 
in bold, unstable ones in dashed-lines. 



plot has a similar aspect and is identical to the left curve 
in figure ||). We will use the notation a x (E) and (12(E) 
to designate branches "I" and "2" (the two roots of 
E 2 = f(a)/ah(a)). 



A. Imposed stress £ 

The evolution of a homogeneous system at fixed E = a 
is truly simple as the Maxwell model yields T = aE which 
can be injected in the equation describing the fluidity, so 
that the evolution of the system is completely described 
by the sole differential equation: 



da dVs . . 
dt da 

where the "effective potential" Vs is : 



V s (a) = / da'[f{a') - a'h(a')E 2 } 
Jo 



(7) 



(8) 



The above equation describes the relaxation of a in the 
potential Vs. Depending on the value of the imposed E, 
this potential has either one or three extrema: 



for E < (J,, 



a 



is the only minimum, 



- for ctm > E > a rn , Vs has two minima, a = and 
a = 02(E) and a single maximum a = ai(E), 

- for E > (Tm , ci = 02(E) is the only minimum. 

The linear stability immediately results (see Fig. ||), 
as the system falls in the potential Vs(a) along the in- 
stantaneous slope at point a(t) (see Fig. g). 

Branch "0" is stable for E < om so that the system 
evolves to this frozen situation whenever the stress is 
too weak E < a rn . At long time the system typically 
creeps but does not flow T(t) = a(i)E with a(t) decreas- 
ing (typically algebraically in time) to 0. For example, 
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FIG. 5: Effective potential Vs(a) for the evolution of a(t) in 
the homogeneous case. Depending on the actual value of E, 
Ve has one or two minima. 



with the explicit values of / and h proposed above, one 
has a(t) ~ 1/t so that T(t) ~ T,Log(t). 

Branch "1" is always unstable, so that for um > E > 
a m the system fluidity evolves either to (again leading 
to creep not steady flow) if the initial value of a is smaller 
than ai(E), or towards a 2 (E) in the opposite case. Note 
that such a behaviour is rather similar to the recently 
reported "viscosity bifurcation" |p2[ . This similarity re- 
mains if the constraint of homogeneity is lifted as will be 
discussed in section IV. The stable branches "0" and "2" 
correspond to stable respectively non-flowing and flowing 
states of the system under E. 

For larger stresses E > <jm if the local flow curve is of 
type IV, branch "0" is unstable and the fluidity always 
evolve towards 02(E). 

Eventually we note that for ojvf > E > er m , the ad- 
dition of noise to the equation describing the evolution 
of a could lead the system to select either branch "0" 
or 'branch "2", depending on which one corresponds to 
the lowest value of Ve (the noise could allow the system 
to jump to the more favourable well, although possibly 
after very long times). Comparing these two values leads 
to introduce a value a* that will prove important in the 
following : for a m < E < a*, = Ve(0) < Vs{a 2 ) so 
that branch "0" is the absolute minimum. For E > a* 
branch "2" is the minimum. Formally, a* is given by 
V a * (a 2 (cr*)) = (see Figure |). 



B. Imposed shear rate f 

Although the steady-state solutions are obviously the 
same as in the previous case, the dynamics is more com- 
plex when r is fixed, as it is governed by the coupled 
equations for the evolution in time of E(i) and a(t). 

The stability of the steady-state solutions can however 
be examined by linearisation around these solutions. The 
two stable branches (0) and (2) on figure [| described in 



the case of an imposed stress remain stable. The stability 
of the decreasing branch (1) at fixed shear rate and with 
homogeneity imposed depends on the choice of the func- 
tions / a nd h in equation (|^) (recall that as mentioned 
in section II B we have formally suppressed the mechan- 
ical inertial mode of destabilisation). For general / and 
h, one can not preclude a destabilisation of the steady 
solution (1) to a limit cycle p3fl. However, this is not the 



case with the expressions of section pi A /(a) 
and h(a) = r'a + a 2 , for which the branch (1 
stable (again when homogeneity is enforced). 



ra 



IV. HETEROGENEOUS FLOW AT IMPOSED 
STRESS E. 

The situation at imposed fixed stress can be appre- 
hended rather simply starting from the analysis of the 
homogeneous case mentioned above with its two stable 
branches "0" and "2" . Heterogeneities are generated by 
the walls which induce wall layers of different fluidities, 
with a major effect on the selection of the steady-state 
flow branch. 



A. Imposed fluidity at the wall 

We start by considering the case where the walls im- 
pose the local value of the fluidity a\ z=0 H — a w . The 
evolution of the system is described by that of a(z, t) 
from which the shear rate can easily be deduced 7(2, t) — 
a(z,t)Ti. Formally, 



da 
~dt 



Sa(z) 



(9) 



which is a rather classical equation describing the local 
relaxation of a "free energy" Fs(a(z)). Here this func- 
tional is: 



H 



1 



[Vv(a) + -D(d z aY]dz (10) 
1 



where it is implicitly understood that a equals a w on the 
boundaries. 

The evolution can be understood as a minimisation 
of Fe, as G?t-Fs = — / dz(d t a) 2 < 0. Schematically, 
minima correspond to adopting throughout the sample a 
value of a minimising Ve (i.e. or 02(E)), and then con- 
necting at the boundaries to the imposed value a w with a 
minimal penalty obtained as a compromise between the 
Ve term and the square gradient term (a more formal 
description of these steady-state solutions is given in the 
Appendix). This leads to steady-state branches in the 
(E, r) plane that deviate from the branches "0" and "2" 
only by the contribution of these wall layers, that shifts 
the average shear rate by a small term proportional to 
i/H, where I ~ {Da w / f{a w )) 1 / 2 is the thickness of these 
layers. 
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FIG. 6: Thresholds values a% and o",j of the stress as a func- 
tion of the fluidity at the walls a w . Above <7; the frozen branch 
is unstable, below ad the fiuidised branch is unstable. When 
the imposed macroscopic stress E is between a d and Oi , both 
branches are stable. Which one is selected depends then on 
the initial conditions, i.e. on the system's history. Areas in 
grey correspond to bulk destabilisation ( "spinodal" ) , whereas 
a simple thick line denotes destabilisation from the wall ( "het- 
erogenous nucleation" ) 
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FIG. 7: Macroscopic flow curve at imposed stress E (thick 
lines). The thin lines (dashed and continuous ) correspond to 
the local flow curve of Figure 2. The value of the fluidity a w 
at the walls fixes the location of the limit of stability o"i and 
Od of the frozen and fluid branches (see Figure 6). When the 
stress is varied monotonously, the shear rate suddenly jumps 
at these values (thick arrows). 



The stability of these branches is now limited by the 
fact that the walls at a w can act as nucleating sites for the 
other phase/branch. For example the "frozen" branch, 
where fluidity decreases from a w on the wall to almost 
zero in the middle, ceases to be stable if E is increased so 
that it becomes interesting (in terms of minimising Fs) 
to increase a beyond a w at the expense of a section of the 
a ~ central part. The threshold roughly corresponds 
to Vz(a w ) — Vs(0) which defines a maximal value of the 
stress (Ji{a w ) for stability on the pasty branch. Beyond 
this value, the a(z) solution "slips" into the potential well 
of branch "2" . Similarly, the fluid branch ceases to be sta- 
ble upon decrease of the stress E when the latter reaches 
the value a d (a w ) defined by V ad (a w ) = V ad (a 2 (a d )). It is 
easy to check that a m < Od < a* , whereas a* < Oi < <tm ■ 
Typical curves for cr,(a m ) and Ud{a w ) are plotted in figure 
I- 

A second process of destabilisation exists. If the wall 
fluidity has extreme values (either low or high) walls may 
no longer be involved in the destabilisation process which 
is then "spinodal" (grey areas in figure 6). 

The conclusion for the macroscopic steady-flow E(r) 
diagram is thus simple. At fixed E there are two stable 
flow curves: 

- for E < (Ji(a w ) a branch close to the axis corresponding 
to a frozen bulk with two layers of finite fluidity close to 
the walls (as schematised e.g. in Figure 2). This branch 
is indistinguishable from the vertical axis for large values 
of the thickness H . 

- for E > (id(a w ), a fiuidised branch very close to "2" 
(identical in the limit H — > oo), with a flowing bulk, and 
slightly less fluid wall layers. 

In the intermediate range ad < E < erj, both branches 



are stable. The branch actually selected depends on the 
system history and preparation (initial conditions) . This 
is very compatible with recent observations of a "bifurca- 
tion" of the asymptotic viscosity at a given value of the 
applied stress for a given preparation 133] . 

If one starts with macroscopic bands of fiuidised and 
frozen material, the fronts between these zones will move 
so that one typically ends on the fluid branch for E > a* , 
and on the frozen branch for E < a*. 

If the stress E is varied slowly in a systematic way (i.e. 
increased or decreased), one has the simple hysteretic 
behaviour described by the thick arrows in Figure |?]. 

The departure from the homogeneous case described 
in the previous section is due to the walls, which intro- 
duce the parameters H and a w . Although the former 
has a moderate role in shifting by a small amount the 
branches in the diagram (by a term proportional to £/ H ), 
the latter controls the position of the instability thresh- 
olds, and thus the amplitude <Ji — ad of the hysteresis 
loop observed macroscopically (actually looking into de- 
tails, these thresholds stresses depend very slightly on H 
but this effect is almost invisible numerically. See Ap- 
pendix) . 

The above picture (Figure is borne out by numerical 
simulations of the equations for various initial conditions, 
and also submitting the system to ramps of increasing 
and decreasing stress. These simulations also allow to 
follow the destabilisation process. Starting from a frozen 
solution, the bulk is destabilised when E = <x;. This 
destabilisation process is initiated at the wall, where the 
fluidity is fixed and non vanishing (acting therefore as 
a destabilising nucleus) : two fronts of large fluidity are 
first created at each wall, then propagate in the material 



8 



2.5 




Z coordinate 

FIG. 8: Destabilisation of the 'frozen bulk' by heterogeneous 
nucleation. The initial state is plotted with dotted lines. As 
the stress is increased so that E > a i} two fronts come from 
the walls, as shown by the arrows before they merge in the 
centre. The steady state is the fluidised bulk. In this simula- 
tion f(a) = a 2 + a 4 , ft (a) = 0.05 *a + a 2 ,a w = 0.8, E = 1.45 



to merge eventually at the middle of the Couette cell to 
give the homogeneous 02(H) bulk solution (Figure ||). 

On the other hand, starting from the fluidised branch, 
when the stress is decreased this highly sheared solution 
is destabilised for S = <Td- The destabilisation process 
depends on a w . If a w < a m , then ad > o m and the 
destabilisation occurs by fronts merging, as previously 
described, whereas in the reverse case ad — cr m , and the 
bulk is destabilised as a whole, as a plateau decreasing 
in time towards zero. 

Note that the problem is formally very similar to that 
of wetting transitions in thermodynamics: it reduces to 
the minimisation of a free energy functional that is the 
sum of squared gradient term (d z a) 2 and a potential with 
two wells Vv(a), with boundary conditions on the walls 
that can favour one or the other of the 'phases'. 

In this interpretation, the sections a* — > o; L on the 
pasty branch, and <Jd - > c* on the flowing one would 
be called metastable. The absence of noise in the equa- 
tions used here makes them absolutely stable. Adding 
some noise in the bulk equation for a would reduce the 
metastable branch, bringing the system (possibly after 
prohibitively long waiting times) to jump from the pasty 
to the fluid branch at er*. 

Note that it may be experimentally difficult to detect 
the hysteresis loop if its amplitude is small (either be- 
cause <Ji — <Jd is small or because of the presence of noise), 
which would make this macroscopic flow curve very simi- 
lar to a textbook description of a fluid with a yield stress 
a*, with possibly a rather abrupt change of the steady- 
state shear rate at this value. 



B. No flux at the wall : d z (a) = 

If the boundary conditions are such that d z a = on 
the walls, then the behaviour of the system is exactly 
that obtained in the analysis of the homogeneous model: 
the system tends to relax any inhomogeneity, and its sta- 
ble steady-state solutions are those of the homogeneous 
model: a pasty branch "0" for £ < om an d a fluid branch 
"2" for £ > cr m . Both branches are stable solutions for 
<j m < £ < om- A hysteretic loop results of amplitude 
larger than in the previous case. 

So the picture is still quite similar to Figure 7 above 
with the following differences: here <7j = um and Od = 
a m , and the branches follow strictly the local flow curve, 
with no slight deviation due to wall layers. 

Again we have checked the validity of this picture nu- 
merically. 

V. HETEROGENEITIES AT CONTROLLED 
SHEAR RATE f 

We now turn to the case where the global shear rate 
is controlled, i.e. the velocity V of the top wall is fixed, 
so that the average shear rate T is imposed. The dy- 
namics of the system is now more complex and governed 
by two coupled integro-differential equation, that can be 
obtained from @ : 

f d t E = -{a) (t)S(t)+r 

\ d t a = -f(a) + ah{a)Y?{t) + ft(o)£<9 t £ + Dd 2 zz a 

(11) 

where (a) is the spatial average of a(z, t), 

(")(t) = -gj dza{z,t), (12) 

and where the boundary conditions are either a| z= o,H = 
a w , or d z a\ z= Q,H = 0. The local shear-rate is given by 
j(z,t)=t+(a(z,t)-{a)(t))Z. 

This set of dynamic equations leads to a much richer 
set of situations, as the dynamics of a does not relax a free 
energy (formally, d t a — —SF^n^ j '5a(z) + ^h{a)dt^ 2 ). As 
a result, for some values of T, the system can have many 
steady-states corresponding to different values of S, and 
for some others none (the system then tends to a limit 
cycle corresponding to a periodic solution for a(z, t) and 
for £(£) that yields a macroscopic stick-slip behaviour). 

To explore this ensemble of scenarios it is necessary to 
solve numerically these equations. Our numerical explo- 
ration leads to a complex picture. Starting with the case 
of a fixed fluidity at the walls, we first describe the main 
features observed, before pointing out more complex fea- 
tures, including the metastability of various patterns of 
shear bands and the occurrence of localised oscillations. 
At this stage we have no analytical understanding of some 
of these features, and of their dependence on the param- 
eters. The behaviours reported here result of numerical 



9 




o o 

y, y m y d 

FIG. 9: Schematic macroscopic flow curve at imposed shear 
rate F (thick lines). Again the thin lines (continuous and 
dashed) denote the local flow curve. The value of the fluidity 
at the walls fixes the location of the limit shear rates 7i and 7d 
of the frozen and fluid homogeneous branches. Between these 
values the system tends to segregate into shear bands. At the 
transition between the two (grey areas), stable or oscillating 
thin layer structures are observed in the vicinity of the walls. 



simulations with / and h fixed to f(a) = ra 2 + a 4 , and 
h(a) = r' a + a 2 . At the end of this section we turn to the 
case of zero fluidity gradient at the walls to underline a 
few differences with the previous case. 



A. Imposed fluidity at the wall. 

We start by analysing the case where the fluidity at 
the walls is fixed to a value a w . Depending on T and a w 
different behaviours are observed. The main features are 
the following: 

- The frozen < X < <Ji and fluidised ad < E branches 



found at fixed stress (section IV), are also found to be 
stable at imposed T. 

- Between two values 7, and 7^ defined in terms of the 
threshold stresses on these branches (see fig. the 
most classical long time behaviour (which we take as 
steady-states) , are solutions corresponding to the coexis- 
tence of "bands" of frozen (a ~ 0) and of fluidised phase 
(a ~ a 2 (a*)) coexisting at a stress E ~ a* . 

- In a narrow vicinity of the transition from the quasi- 
homogeneous branches to this situation of coexisting 
"shear bands" (i.e. T slightly larger than 7, or slightly 
smaller than 7^), various phenomena can occur in the 
vicinity of the wall: nucleation of thin "layers" of fluid- 
ity different from the bulk, or of "layers" of periodically 
oscillating fluidity. 

These behaviours have different signatures as shown 
on the macroscopic flow curve sketched on figure |9|. 



1 . Frozen and fluidised regimes 

First, for small and large imposed shear rate, T one 
recovers the solutions obtained at imposed stress : except 
for small boundary effects, they correspond to the frozen 
branch " 0" , and to the fluidised branch "2" of the local 
flow curve. 

Along the fluidised branch, the E(F) relationship al- 
most exactly follows the intrinsic flow curve obtained in 
the homogeneous case. The frozen branch actually dis- 
plays an almost vertical section going to zero in the limit 
of vanishing shear rate T — » due to the presence of a 
small wall layer of finite fluidity imposed by the wall, a w , 
and of thickness £. 

These domains are limited by the global stresses ai(a w ) 
and <Jd{a w ), as discussed in section [V. The correspond- 
ing shear rates 7^ and jd can be calculated numerically 
with the method described in the Appendix. 



2. Shear-banding regime. 

A shear-banding domain connects the two previous 
frozen and fluidised regions. It corresponds to a " coexis- 
tence" in the shear cell of the frozen and of the fluidised 
states. Its signature on the macroscopic flow curve is 
a quasi-plateau at a stress E ~ a*. All over the band- 
ing domain (for a given a w ), the stress is found to re- 
main in the vicinity of the value a* defined in section III 
by Vcr* (0,2(0*)) = V a *(0). This approximately fixes the 
value of the fluidity in the fluidised region to a* = ai (cr*), 
with 02(17) being the correspondin g br anch obtained in 
the homogeneous case (see section QI). The shear rate 
in the fluidised region is approximately given accordingly 
as 7* = a* a* . The total width of the frozen and fluidised 
regions are determined by the global constraint on the 
total shear rate. The width of the fluidised region hf 
is typically given by hf ~ HT/j*, with 7* defined here 
above (small boundary effects neglected). 

An interesting question concerns the number of flu- 
idised bands : how many fronts between a frozen and a 
fluidised bands can be stabilised for given H, a w , and 
F ? In order to answer this question, a 'zoom' of the 
quasi plateau is plotted on figure |l0| (a) for given H and 
a w . The stable quasi plateau actually consists of parts 
of the branches representing the steady states solutions 
with several fronts (on the figure these stable parts of 
the branches are plotted with thick lines). The steady 
state multi front branches were calculated numerically 
with the method detailed in the appendix. The fluid- 
ity profiles corresponding to the branches with one and 
two fronts are represented on figure [lO] (b). Their sta- 
bility was then examined with the numerical simulation 
described in II D. 



No analytical characterisation was found to the limits of 
stability of the multi front branches. Beyond these limits, 
either oscillations are observed, or other stable branches 
are reached. Their complex dependence on a w and E, 
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FIG. 10: (a) Calculated steady state flow curves : the unsta- 
ble stationary solutions are plotted with thin lines, the stable 
solutions with thick lines, the oscillations with dotted lines. 
The numbers on the branches represent the number of fronts 
between the fluidised and frozen bands, r = 1, r' = 0.01, H = 
40, a w — 0.25. b) Fluidity profiles corresponding to the points 
marked with crosses on figure (a) for the branches 2a,2b, and 
3. 



is shown on figure O. Note that, with the parameters 
taken to plot figure |10| (a), for a small range of T (just 
left of the middle cross in Figure 10a), two solutions are 
possible (depending on the history of the system) : sta- 
ble oscillations of a central fluidised layer (branch 2b), 
and a steady solution with a single fluid layer at the wall 
(branch 1). In the case of oscillations, when T is further 
decreased to the limit of stability of branch 2b, the sys- 
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Shear rate 



FIG. 11: Behaviour at low shear rate for a fixed H — 100, 
r = 1, r' = 0.01. The solid line corresponds to a w i — 0.3 and 
the dashed line to a W 2 = 1. In the case of a ml , oscillations 
are observed, and the envelope of the oscillations amplitude 
is represented 



tern evolves toward the pattern of branch 1 as indicated 
by the arrow in Figure 10a. 



3. Stick slip behaviour 

When no stable branch is reached, a time dependent 
oscillating behaviour is evidenced for stress and fluidity. 
As seen on figure occurs for small wall fluidity and 

imposed shear rate. This behaviour is very reminiscent 
of stick slip, a phenomenon commonly encountered in 
solid friction []33[ |34f . Here it originates in an oscillating 
behaviour of the fluidity between a " frozen" state (" stick" 
phase) and a partially fluidised state ("slip" phase). 

The fluidisation however occurs only within one 
layer, the rest of the system being frozen. This spa- 
tially localised oscillating behaviour resembles a periodic 
fracture-healing process, as observed e.g. by Pignon et 
al. jn^. The fluidised layer is localised either close to 
the wall, or in the middle of the shear cell, its fluid- 
ity oscillates between two extremal values, and its thick- 
ness is approximately fixed. The time evolution of S and 
the maximum fluidity of the layer are plotted on figure 
|l2] , showing characteristic features of a "stick-slip" be- 
haviour. 

The system rapidly relaxes to the frozen state, i.e. a 
very small fluidity, while the stress increases until a limit 
where the fluidity abruptly increases, i.e. the fluid is 
suddenly sheared. This sudden increase is due to the 
penetration of the maximum fluidity in the well potential 
of ci2{<j), that is enlarged as a increases. To maintain a 
fixed global shear rate, a relaxes quickly and the cycle 
starts again. 
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neous fluidity in the cell. Moreover, as discussed in sec- 
tion III, the decreasing branch is linearly stable at fixed 
shear rate (remember that we have suppressed inertial 
terms), and no destabilising mechanism is induced here 
by the boundaries (in contrast with the fixed wall fluid- 
ity case). Formally, the system is thus attracted either 
by the homogeneous solution or by the inhomogeneous 
one, depending on initial conditions (an initial inhomo- 
geneous state often leading to a banded inhomogeneous 
solution) . 

The macroscopic flow curve for (11,12) with the 
present boundary conditions consists of: 

- the two increasing branches of the local curve ( the 
branches are parts of the local curve as the solutions 
are homogeneous). Those branches are also stable 
in the case of imposed stress. 

- the decreasing branch of the local curve. This is 
of course an artefact due to our neglecting iner- 
tial terms, which suppresses the "mechanical insta- 
bility" that always destabilises this branch. 

- the quasi plateau, with £ ~ a*. As in the case of 
an imposed fluidity at the wall, it is constituted of 
parts of multi front branches, that become unstable 
at different values of the shear rate. 

- an oscillating domain, for small shear rates, when 
the single front branch looses its stability. 

Ruling out the unphysical decreasing branch, the picture 
is thus quite similar to Figure 9 (without the effective 
wall slip in branch 1 due to the finite wall fluidity there). 
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FIG. 12: The amplitude of the oscillations of the global stress 
and of the fluidity of the oscillating layer are represented as 
a function of time (a), and in the plane (a, E) (b), in the case 
a w = 0.3, t = 1, r' = 0.01, f = 0.12. Its width hardly varies. 



B. No flux boundary condition at the wall 



o 



In the case of a no flux boundary condition and for 
intermediate values of Y (7, < Y < 7^), two very differ- 
ent locally stable sets of solutions of (11,12) are possi- 
ble : one is very similar to the solutions obtained above 
for fixed fluidity at the walls and exhibits shear banding 
and "stick-slip" behaviour, the other one is the homoge- 
neous solution at fixed shear rate obtained in section 111 . 
This peculiarity is associated with the specific form of the 
boundary condition, which is compatible with a homoge- 



C. Comment on oscillations 

Oscillating response under steady driving has re- 
cently been numerically observed in a model for shear- 
thickening materials |3q] , which in contrast to the present 
one described the homogeneous behaviour of a system (no 
spatial heterogeneity in the formalism), and the oscilla- 
tions were shear-rate oscillations at fixed stress. Formally 
there are some similarities as in both cases a continuous 
(and thus infinite) set of variables evolves through cou- 
pled differential equations with non-monotonic functions 
involved, and driven by a global integral constraint. How- 
ever, in contrast to the situation of ]3j| where all the vari- 
ables were involved in the oscillation, here only a limited 
set of variables do (corresponding to a spatially localised 
oscillating layer). 

In principle a finite number of variables is enough to 
get an oscillating rheological response (limit cycles) under 
steady driving [|36| , so it is no surprise oscillations show 
up here given the complexity of the set equations ( 1 1 p"^ ) 
and the number of variables involved. One may even ex- 
pect more complex dynamical behaviour such as some 
form of chaotic response |36| , but we found none within 
the set of parameters numerically explored. Clearly, a 
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truly three dimensional version of this model, with the 
corresponding non- linear convective terms, should in- 
crease the likeliness of turbulent/chaotic behaviour (even 
in the absence of inertia |B7fl). 



VI. DISCUSSION 

In the present paper we used an extension of a model 
previously introduced to account for the rheology and ag- 
ing in pasty systems, in which the steady state local flow 
curve is non-monotonous so that heterogeneous flows can 
be generated at steady state. The key ingredient of the 
model is a coupling between the mechanical properties, 
namely shear stress or shear rate, and local relaxation 
within the fluid, characterised by a "fluidity". Although 
analytically very simple, the model yields a rather rich 
macroscopic phenomenology. 

• Under controlled stress E, the steady state be- 
haviour is schematised on figure The fluid un- 
dergoes an hysteretic phenomenon if the stress is 
increased or decreased progressively. The effective 
yield stress is shown to depend on the fluid history 
and on the wall conditions. Besides, no banding is 
observed and there is a discontinuity in shear rate 
r which occurs at specific values of the stress. 

In this case, we have shown that the dynamics re- 
duces to the minimisation of a free energy function. 
This is generalisable to any function /, g describing 
aging and fluidisation, and to more complex diffu- 
sion terms D(a). Note that the dynamics was not 
built as a free energy relaxation but naturally ap- 
pears this way (see related discussions in |ll|). This 
gives some generality to the qualitative results, be- 
yond the specific model studied here. 

As no noise term was included, the equivalent of 
homogeneous nucleation in thermodynamics prob- 
lem was suppressed and boundary conditions play 
a crucial role in destabilisation processes. Hence, 
the values of the thresholds depend on the bound- 
ary conditions. However, the qualitative behaviour 
(hysteretic loop) is rather insensitive to the specific 
type of boundary condition chosen. 

• At fixed global shear rate T, the macroscopic be- 
haviour is schematised on figure : banding solu- 
tions show up, separating a "frozen" region from 
a fluidised region. The shear rate within the fiu- 
idised band is hardly dependent on the imposed 
global shear rate, and fixed only by characteristics 
of the material. Within the banding domain, the 
thickness of the fluidised region increases (nearly) 
linearly with the global shear rate T to eventually 
fill the full shear cell. 

At small (total) shear rate, an oscillating behaviour 
of the stress E is observed, with features very rem- 
iniscent of a classical "stick-slip" behaviour. The 



system is globally in a frozen state, and fluidises 
periodically within a thin layer. Such a behaviour 
compares well with the fracturing/ healing pro- 
cess observed by Chen et al. || and Pignon et 
al. Gd|. These localized oscillations are to be con- 
trasted with homogeneous oscillations obtained in 
other models for different systems (e.g. pq]). 

Our observations result from numerical simulations 
and the influence of the specificity of the modeli- 
sation of the fluidisation and aging remains to be 
understood. Again, the different thresholds de- 
pend on the boundary conditions, but the quali- 
tative behaviour described are robust to changes in 
the types of boundary conditions. Extension of the 
present model to a fully three-dimensional descrip- 
tion could however yield a more complex dynamical 
behaviour. 

In conclusion, we suggest that analysing yield stress 
fluids close to the yield stress at both constant stress and 
constant shear rate could provide interesting information. 
It could also prove informative to check the sensitivity of 
the results to the nature of the surfaces bounding the 
fluid. 



APPENDIX A: APPENDIX : STATIONARY 
SOLUTIONS OF EQUATION @ 

In this appendix, we detail the steady state solutions 
of Eq. (Rjl) recalled here : 



a(z)T, = j(z) 

-f(a(z)) + a{z)h{a{z))Y? + Dd 2 zz a{z) = 
a\ z =o,H = a w or d z a\ z = ,H = 



(Al) 



For clarity, the potential Ws(a) is introduced as : 

We(o) = -J a da>[f(a') - a'h{a')Y?]. 

Note that Wy, = — Vk as defined in Eq. (||). Hence for 
a given E, we look for a solution a(z) such that 



-d a Wz(a) = Dd 2 zz a(z) 
a\ z =o,H = a w or d z a\ z=0 ^ H 







(A2) 



the local shear rate is then immediately deduced 

7(2) = a(z)£ 



Equation (A2) can be interpreted as the classical equa- 
tion of motion of a particle of mass D and position a in 
the potential Wy. where z plays the role of the time. The 
energy E m = \{d z a) 2 + Ws{a) is conserved during the 
motion. 

With this interpretation, the boundary conditions im- 
pose either initial and final positions (a w imposed at the 
walls), or velocities (d z a = imposed at the walls). 

We choose to detail here the case of an imposed a w at 
the wall, and the stress E is chosen <7j > E > er^ (see 



section IV) 
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FIG. 13: a) Potential W s . E = 1.44, r = l,r' = 0.01. b) 
Fluidity profile for a fluidised bulk. a w — 1, H — 15. 



The potential is represented on figure [13|. The ex- 
treme values of the 'mechanical energy' E m are Ws(a w ) 
and Wz (0,2)- From the analysis of these solutions, it is 
possible to progress further in answering a few questions 



The fluidity profile is represented on figure [13| Note 
that E m = We (02) corresponds to a fluidised bulk, 
with an infinite thickness H ~ 00. 

With d z a\ z= o > 0, a profile with a bulk of fluidity 
higher than a w is described. 

Symmetrically, a profile on the frozen branch, that 
is with a bulk fluidity smaller than a w is described 
by fixing d z a\ z= o < 0, and following the same pro- 
cedure. 



• What is the limit of existence for the , 
solution Oi ? 



luidised bulk 



ai(H,a w ) was defined in section IV as the stress 
beyond which a profile with a frozen bulk can not 
be observed. 

For H ~ 00, ai(co,a w ) is defined by W ai {a w ) — 
0. With the mechanical analogy, the mechanical 
energy E m is set to 0, that is the particle leaves 
its position a w , with a velocity d z a\ z =o ~ 0, and 
spends an infinite time at the position a m = 0. 

For a > (Tj, W a (a w ) > 0, and for a < a w ,W rr (a) < 
W a (a w ). With the mechanical analogy, the parti- 
cle can not leave the position a w with a negative 
velocity and come back to it after a time H . There 
is no solution with a frozen bulk. 

For a finite H, a frozen bulk requires E m < 0, there- 
fore o~i{H, a w ) is slightly smaller than Cj(oo, a w ), by 
an amount that scales as H~ 3 . 

The situation is symmetric for the definition of 
o~d(H,a w ), which is found slightly larger than 
cr d (oo, a w ). 

• How can one calculate the fluidity profile for a so- 
lution with several bands? 

The mechanical energy E m is now decreased from 
Mm(W E (0), Wx{a 2 )), to W^{a w ) to find banded 
stationary solutions : a frozen bulk solution with 
thickness h\ and, a fluidised bulk solution with 
thickness h 2 (both with the same energy E m ) such 
that hi + h 2 — H form a single interface solution. 
The same process is derived to find multibanded 
solutions. Solutions of figure |l^ were calculated 
numerically with this method. 



• What is the fluidity profile on the frozen and flu- 
idised branches of figures Q andfy? 

For a given thickness H, E m is set so that a| z=0 ,ff = 
a w with d z a\ z=0 > 0. With the mechanical analogy, 
this amounts to fix the mechanical energy or initial 
velocity for the particle to leave position a w , reach 
a maximum a m with d z a\a m — 0, and go back to 
its initial position within a time H. 
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